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Abstract 

m 

J> ' We analyze lepton flavor violating r — > ^7 and /i — > e7 processes in SUSY 

GUT model in which sfermions have special mass spectrum. It is assumed that 
only third generation sfermions which are contained in 10(Q,U C ,E°) of SU(5) can 
have a different mass from the others. This mass spectrum is led from Eg GUT 
model with horizontal symmetries. It is shown that branching ratios of r — > /X7 
and \i — ► e7 depend strongly on a right-handed stau mass. The weak scale stability 
requires the light stau, so large decay rates can be expected in this scenario. When 
the stau mass is around 150 GeV and tan (3 ~ 10, the branching ratios can be larger 
than Br{r — > ^7) ~ 10 -8 and Br(u — > e7) — 5 x 10 -12 , which are within reach 
of future experiments. In addition, this model has an interesting feature that the 
final state charged lepton tends to have the right-handed chirality. 



1 Introduction 

The measurement of the Lepton Flavor Violation (LFV) in charged lepton sector is one of 
the most important measurements to search for New Physics beyond the Standard Model 
(SM). Since the neutrino oscillations have been observed, the charged LFV processes can 
be expected to be non-vanishing. However, the SM with the Majorana neutrino masses 
leads too small LFV decay rates to be observed due to smallness of the neutrino massespQ. 
Thus, detections of the LFV processes means discovery of New Physics. 
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In fact the supersymmetric extension of the SM leads generically large LFV decay 
rates, because of picking up the off-diagonal entries of slepton mass matrices (in substitu- 
tion for neutrinos) in intermediate state. However, it is well known that Flavor Changing 
Neutral Current (FCNC) and CP violating processes are highly suppressed experimen- 
tally. They impose severe constraints on the minimal supersymmetric (SUSY) standard 
model (MSSM) parameter space[2j. The simplest way to satisfy these FCNC constraints 
is to adopt the universal sfermion soft mass matrices at some scale. 

In many analysis that have been done so far, the universal soft masses at the gravi- 
tational scale or at the grand unified scale have been assumed and the LFV decay rates 
are calculated by the non-universality which is induced by the radiative correction^ Hj. 
For example, in the MSSM with the Majorana neutrino masses, the radiative correction 
induces the non-universal sfermion masses [3]. Unfortunately, it is difficult to predict the 
LFV decay rates quantitatively in this scenario unless the Dirac neutrino Yukawa cou- 
plings are fixed, because such LFV decay rates strongly depend on the Dirac neutrino 
Yukawa couplings, which looks impossible to be measured directly. In SO (10) grand uni- 
fied theory (GUT), the Dirac neutrino Yukawa couplings are related with the up- type 
quark Yukawa couplings. Then we can predict the LFV decay rates [3], but they are too 
large to be consistent with the experimental bound if the diagonalizing matrix of the 
fermions mass matrices for doublet lepton has the Maki-Nakagawa-Sakata[5j(MNS)-Me 
large mixings. In this direction, some special structure which suppresses the // — > 
process is required. 

The assumption of universal soft masses is not necessarily required. For example, 
since the FCNC constraints for the first two generation fields are much more severe than 
for the third generation fields, the universal soft masses sometimes have been imposed 
only for the first two generation sfermion, which are realized if we impose non-Abelian 
horizontal symmetry, for example, U(2), under which the first two generation fields are 
doublets and the third generation fields are singlets[6]. However, if the diagonalizing 
matrix of the fermions have the MNS- like large mixings, such sfermion mass spectrum 
leads to too large FCNC to be consistent with the experimental bound. Therefore, such 
non-universality must be introduced only for the sfermion whose fermionic superpartner 
has the diagonalizing matrix with the Cabibbo-Kobayashi-Maskawa[7](CKM)-/z&e small 
mixings. 

In the SU(5) SUSY-GUT, it is a reasonable assumption that the diagonalizing matrices 
of 10 representation fields have the CKM-like small mixings, while those of 5 fields have 
the MNS-like large mixings, because 10 representation includes the doublet quarks, and 5 
includes the doublet leptons. Under this assumption, only the third generation sfermions 
of the 10 representation can have different masses from the others without contradicting 
with experimental constraints from various FCNC processes. In this paper, we study the 
LFV under the non-universal sfermion masses above. 

Before examining the LFV, we summarize several characteristic features of the non- 
universal sfermion masses. First of all, the rotation matrix for 10 which makes fermions 
to mass eigenstate is expected to have the CKM-like small mixings, since 10 involves 
quark doublet. Therefore, large off-diagonal entries for 10 representation sfermion mass 
matrix do not arise after this rotation, even if the initial soft masses are not degenerate. 
Moreover, FCNC constraints among 1-3 or 2-3 generations are not so severe than that 
of 1-2 generations. Therefore, we can expect that non-degeneracy for third generation 
does not conflict the FCNC constraints. And it is important that we can expect larger 
FCNC for the third generation fields than in the usual universal sfermion mass case. 
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Second, the naturalness of the Higgs mass in MSSM requires that the gaugino masses, 
the higgsino mass, and the stop masses must be around the weak scale. In the above non- 
universality case, we can take the larger masses for the other sfermion masses than the 
third generation sfermions of the 10 representation without conflicting the naturalness 
arguments, because the both stops (left and right handed stops) are included in the third 
generation 10 representation of SU(5). In principle, we can take such mass larger than 

1 TeV, which can weaken the various constraints from FCNC, electric dipole moments 
(EDM) and the muon g — 2, etc. related with the first two generation fields. Again, the 
FCNC related with the third generation fields can become relatively large, which may be 
detectable in future experiments. Finally, such non-universality of the sfermion masses 
is naturally realized in the E 6 SUSY-GUT modelJHl IS] with SU(2) or SU(3) horizontal 
symmetry [in]. 

In this paper, we will analyze the lepton flavor violating r — > /ry and \i — > processes 
under this non-universal sfermion masses. The point is that the naturalness requires that 
the 10 3 sfermion masses are 0(A W ), i.e., the right-handed stau mass becomes 0(A W ). 
We show that if mass of right-handed stau is near the weak scale, each LFV decay rate 
becomes sufficiently large to be found by near future experiments. Actually, the LFV 
decay rates are very sensitive to the right-handed stau mass. Thus, if r — > ^7 or fi — > cy 
processes are discovered and the branching ratios are measured, this model may predict 
a mass of right-handed stau. We will see that this model has interesting feature, that 
is, this model indicates that the final state charged lepton tends to be right-handed, 
since main source of flavor violation comes from the Feynman diagram which includes 
the right-handed charged slepton. 

This paper is organized as follows. In the next section, we will specify the model, 
and discuss how lepton flavor violating processes take place in this model. The results 
of our numerical study are given in Section 3. In Section 4, we will give a summary and 
discussion of our results. 

2 Lepton Flavor Violation in this model 

In this section, we discuss how lepton flavor violating decays take place in a SUSY-GUT 
model with the non- universal sfermion masses mentioned in the previous section. First, 
we assume that the SM gauge group SU(3)c X SU(2)l x U(1)y is unified into SU(5). 
(Note that the following arguments can be applied to unified models with larger unified 
gauge groups, for instance, SO(10) or E 6 .) Then, around the GUT scale, the Lagrangian 
for the SUSY breaking terms is given as 

— C-soft breaking = —£-soft mass + (^ijlOilOj^H + ^ 6 10j5j5# 

+ AZjUtfiH + B5 H 5 H + ^M 1/2 X a X a + h.c.) , (1) 
-Csoft mass = (m^yioj Wj + (m^ijtfij + ml^jUlj + m 2 Hd 5 i H 5 H + m 2 Hu & H 5$) 

where the tilded fields are regarded as sfermions, A a is SU(5) gauginos, 10j = (Qi, Ui, Ei), 
5i = (Di,Li), h = Ni (i = 1,2,3), M x p is SU(5) gaugino mass, and -C soft mass gives 
sfermions and Higgs soft mass terms. 5# and 5# include the down-type Higgs H d and 
the up- type Higgs H u , respectively. 

As mentioned in Introduction, we consider following forms for the sfermion mass 
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matrices. 




The reasons are as follows: 



• Severe Constraints on the first two Generation: 

There are severe constraints on the off-diagonal entries of low energy sfermion mass 
matrices from the various FCNC processes in the basis of the superfields in which all 
quark and lepton mass matrices are diagonal ( super-CKM basis). This is famous as 
the SUSY-Flavor problem. In particular, the (1,2) component suffers from strong 
constraints by CP violating parameters of K°-K° mixing, lepton flavor violating 
(i — > ej, and so on. Thus, we cannot release degeneracy between 1-2 generation, 
otherwise a large (1,2) element arises in super-CKM basis and it conflicts with 
experiments. 

• Large Neutrino Mixing and Constraints on 2-3 Generation: 

It is known that there is a maximal mixing between 2-3 generation in lepton 
sector [TTJ. This large mixing sometimes is related with the mixing of 5 since 5 
contains left-handed lepton doublet. Thus, if we release degeneracy between 2-3 
generations in 5, a large 2-3 component of 5 is generated in super-CKM basis. It 
conflicts with constraints from B°-B° mixing and r — > fij. Therefore, we should 
assume degeneracy of 2-3 generation for 5 soft masses. While, we can take a soft 
mass of the third generation of 10 different from that of the other generation, since 
a mixing of 10 is expected to be as small as CKM mixing. 

• Naturalness of the MSSM: 

In the MSSM, there is no terrible fine tuning problem due to supersymmetric can- 
cellation for quadratic divergence for radiative correction to Higgs mass square. 
Instead of the quadratic divergence, the radiative corrections are proportional to 
the SUSY breaking mass square, and therefore SUSY breaking mass scale must be 
around the weak scale Ajy to avoid some tuning if the couplings are 0(1). Among 
the various sfermion masses, only the stop masses must be the weak scale, while 
the other sfermion masses can be taken much larger than the weak scale because 
of the smallness of the Yukawa couplings. One of the advantages of the mass spec- 
trum we assumed is that we can take all the sfermion masses much larger than 
the weak scale except the masses of sfermions which are contained in IO3. In such 
parameter region, several constraints in SUSY models, for example g — 2, EDM, 
FCNC, etc, are relaxed without loss of naturalness. Especially, A-term contribution 
to EDM or FCNC for the first two generations can be decoupled in this region. In 
the literature, it is argued that the stop masses are required to be larger than the 
weak scale to realize the Higgs mass larger than the lower bound > 114.4 GeV 
obtained at LEP II. Such large stop masses lead to some tuning in the Higgs sector. 
This is called "the little hierarchy problem" [12] . However, as discussed in [13], the 
stop masses around 300-400 GeV can be consistent with the LEP II experiments, 
if the lightest CP-even Higgs has a small ZZh coupling. Therefore, some of the 
parameter region considered in this paper may lead to such small stop masses, but 
it is consistent with the LEP II results of the Higgs search experiment. 
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• E 6 GUT with U(2) horizontal Symmetry: 

Introducing a non-Abelian horizontal symmetry is one of various interesting pos- 
sibilities to realize the flavor blind nature of the sfermion masses [HJ Q3]. Most 
of the previous attempts adopt U(2) or SU(2)x U(l) symmetry as the horizon- 
tal symmetry [H], which realize the universality of the first two generation sfermion 
masses and the large top Yukawa coupling, if we assign doublets for the first two 
generation fields and singlets for the third generation fields and the Higgs fields. 
Unfortunately, the universality of the first two generation fields is not sufficient to 
suppress the FCNC constraints especially when the diagonalizing matrix of fermion 
mass matrices have the MNS-like large mixings as we pointed out in the previous 
section. E 6 GUT x U(2) Horizontal Symmetry model[TU] succeeds to give not only 
the observing large hierarchy of fermion masses but also the sufficient universality 
of the sfermion masses to avoid FCNC or EDM constraints in a simple way even 
when the diagonalizing matrix has the MNS-like large mixings. It is essential that 
this model can naturally realize the sfermion mass spectrum which we discussed in 
this paper because of an E 6 nature. 

The superpotential is given at the GUT scale as follows. 

W = V,'< 10,10/^ + Y$ e IO^h + Y^l^hjj + ^Mijlilj + fi5 H 5 H . (4) 

Since 10 contains left-handed quark doublet, the diagonalizing matrices of 10 fields are 
naively expected as 





A 


A 3 




1 


A 2 


A 3 


A' 2 





10^(^0)^10,-, Vio ~ A 1 A 2 , (5) 



where A is Cabibbo angle (A = 0.22) and we omit the 0(1) coefficients. This rotation 
generates off-diagonal entries of the sfermion masses of 10 fields as 



2 \ / 2 

m \ / m o 

2 It/., I ™2 | , _2 



Vi T o ™o \ v w ~ m + m 



™30y 







A 5 




A 4 




A 2 





X b A 5 

<«o~K-^o) h 5 A 4 A 2 |, (6) 

\A 3 A 2 

at the GUT scale. 

In the following, we examine the LFV processes with the above sfermion masses. We 
use the mass insertion method to estimate the amplitudes for a while. As the slepton 
mass matrices at the SUSY breaking scale, we adopt the slepton mass matrices given 
by replacing the m , m 30 , and m mix0 with m, m 3 , and m mix , respectively, which are the 
slepton mass parameters at the SUSY breaking scale. As mensioned in Introduction, 
we neglect the radiatively induced off-diagonal slepton mass matrices. The weak scale 
stability requires that m 3 , Mi/ 2 , fi, and A must be around the weak scale, but the 
other parameters, for example, m, can be much larger than the weak scale. In this 
paper, we often take m 3 , Mi/%, fi, A <C m. For the left-right mixing SUSY breaking 
parameters, A, we just neglect the contributions in this paper. One reason is that this 
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contributions are decoupled in the limit 772,3, Ml/2, /i, A <^ m, and the other reason is that 
we would like to examine the model independent parts while the left-right mixing SUSY 
breaking parameters are strongly dependent on the explicit models. On the other hand, 
the diagrams (a) and (b) are not decoupled because m 2 ^. becomes large in the limit. 
Therefore, we consider only (a) and (b) for a while. When tan/? is large, the diagram 




(b) which is proportional to Y ej = \f2m ej jv cos (3 is significant, where v = 246 GeV. The 
ratio of (b) and (a) is roughly 

Amp {b) 

- —~t a n p. (7) 

Amp (a) v ' 

Let us roughly estimate the branching ratio of [i — > e 7 from the diagram (b) which is 
expected to give the largest contribution. Comparing the diagram for /i — > e 7 with that 
for fi — > evv (Br(fi — > evv) ~ 100%), we can roughly estimate the fi — > e 7 branching 
ratio as 

Br(/i^ej) ~ Br(fi — > evv)— — 3- ~ -— j- e(g ) — 



'T(fi-^evv) l6ir 2 G 2 F 
47r«(a') 2 (tan/?) 2 A 10 



m s 



G\ m% 



where Gp is Fermi constant and ms is typical superparticle mass. For example, if we 
take ms = 400GeV, tan/3 = 5, (jHJ) is found as 

Br{fi -> e 7 ) ~ 10" 11 (xO(l) 4 ). (9) 

We find that it is as large as the experimental upper bound (Br(fi — > e 7 ) < 1.2 x 
10~ n [15j). Since there is no lepton mass dependence in equation ([H]), it can be applied 
for the r — > /i 7 branching ratio. Substituting (777^^)32 and Br{r — > evv) for {m 2 mix )2i 
and Br({i — > e^z/), respectively, the r — > yU 7 branching ratio is found roughly as 

„ . . (Brir^evv) A 4 \ _ . . 

~ 10 3 x5r(/i^e 7 ) (xO(l) 4 ). (10) 

Then, for the same ms and tan/? as in the above estimation of Br(fi —>■ e 7 ), 

Br(r -> /x 7 ) ~ 10~ 8 (xO(l) 4 ). (11) 

The branching ratio of r — > /i 7 is also as large as reach of experimental bound [Br{r — > 
/i 7 ) < 4.5 x 10~ 8 P2)]). The relation between Brir — > /i 7 ) and Br(r — >• e 7 ), 

Br{r -> e 7 ) ~ A 2 5r(r -> /j 7 ) (xO(l) 4 ), (12) 
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is also an interesting prediction of this scenario. Unfortunately r — > &y process is more 
difficult to be observed because of the small branching ratio because the experimental 
bound, Br{r — > &y) < 1.1 x 10~ 7 [16j, is similar to that of r — > ^7 process. Since the 
above estimations are very rough, more accurate calculations are necessary. Such accurate 
calculations and the results are given in the next section. 

Before going to the numerical calculations, we mention about the other source of 
LFV. It has been discussed that the loop corrections induce the off-diagonal elements of 
left-handed slepton soft masses [3] as 

Hh - -^( 2 ^o + <J(^")« log , (13) 

where M and M me( nated are the right-handed neutrino mass scale and the mediation scale 
of SUSY breaking, respectively. If 

(Y^Y V ) 23 > X 2 or (Y^Y U ) 12 > A 5 , (14) 

then this loop correction may become sizable. For example, in SO(10) GUT model 
with minimal Higgs sector above equality is just hold because of SO(10) GUT relation 
yv _ yu ^ though not in E 6 models due to the smaller Dirac neutrino Yukawa couplings. 
Since these contributions are decoupled again in the limit m 3 , Mi/ 2 , fi, A -C rn, we 
neglect this effect in this paper. 

In case where A d,e have non- vanishing off-diagonal elements in the Y d,e diagonal basis, 
it leads large LFV. The LFV occurs through the diagram (c). The ratio of amplitudes of 




diagram (c) and that of (b) is roughly found as 

Amp(c) M 1 A d fv cos (3 

~ (15) 

Amp (b) m h tan (3 (m^ji ' 

If A d f > Yf tan/3(m^ lia ,)jj/Mi is hold, a contribution from diagram(c) is more significant 
than (b). As we noted, this contribution is dependent on the explicit models and decou- 
pled in the limit 777,3, M1/2, fi, A <C rn, and therefore, we just neglect this effect in this 
paper. When the A-terms are proportional to relevant Yukawa matrices (A d,e = aoY d,e ) 
or even when the A-term is vanishing, the double-mass insertion diagram (d) via [i param- 
eter can contribute the LFV processes. This contribution is included in our calculation 
in the next section, but it is small. This may be because it is decoupled in the above 
limit and it is vanishing when 7773 = m. 
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3 The Branching Ratios 



Let us discuss the branching ratio of r — > /i7 and /i — > e^y. The amplitude of lj — > If) is 
written as 

T = e"(g)^(p - q)i ( y a ^{Af ) P L + Af P R )u 3 (p) , (16) 

where p and q represent the momenta of lj and photon, respectively. Ignoring a mass of 
li the decay rate is given by 

TYl 

nh^kl) = ^{\Af ) \ 2 + \Af ) \ 2 ). (17) 

Now we can see an interesting feature in this model. This model says final state lepton 
tends to have the right-handed chirality. From the operator form ffTBj) . we find that 
chirality flip is necessary with picking up the Yukawa couplings. In this model, only 
flavor violating coupling in the charged lepton sector is rv? mix . So, the intermediate state 
must have the right-handed chirality. If the final state lepton has the left-handed chirality, 
we must pick up the Yukawa couplings at the final vertex (Yj). On the other hand, if 
the final state has the right-handed chirality (then the initial state has the left-handed 
chirality) we must pick up that at the initial vertex (Yj). It means Al is as large as A R 
times mj/rrii, so we find the final state lepton tends to have the right-handed chirality. 
We can test this prediction experimentally by measuring the angular distribution of final 
lepton. 

In the basis where the charged slepton and the neutralino are mass eigenstates, the 
amplitude of lj — > Ifi is given by 




where mg k (k = 1, 2, • • ■ ,6) and M^o (a = 1, • ■ ■ , 4) denote the masses of the charged 
slepton and the neutralino masses, respectively. Gq c r and H QeR are the effective couplings 
which originate with the gauge interaction and Yukawa interaction, respectively. The 
explicit form of G 0eR and H 0eR and the function F 2 (x), F±(x) are defined in Appendix. 
We assume the GUT relation for gaugino masses 

Mi _ M 2 _ M 3 

2~ 2~ 2~" \^) 

9i 9i 91 

Moreover, we use the observed CKM matrix as the form of unitary matrix for e R which 
diagonalizes Y e 0. Here, only for simplicity, we take the vanishing CP phase. In the 
above assumption, the number of free parameters which determine the branching ratios 
is 5 : m, m.3, the higgsino mixing parameter (//), the SU(2)l gaugino mass (M2), and 
the Higgs VEV ratio (tan/3). In the following, instead of m and m 3 , we use the first two 
generation right-handed slepton masses defined as m| B = m 2 + m| cos 2 2/3 sin 2 6w and 

the right-handed stau mass defined as ?^e H3 = ^3 + rn? z cos 2 2/3 sin 2 9w, which give the 

1 Note that the effect of 0(1) coefficients can be large as discussed in the previous section. We also 
neglect a suppression factor (~ 15%) induced by electromagnetic corrections [17]. 
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eigenvalues of the slepton mass matrix in the limit in which the lepton Yukawa couplings 
are vanishing. 

In Fig. 1, we show the dependence of Br(r —>■ /ry), Br(fi —>■ erf) on the first two gen- 
eration right-handed slepton mass. Here, we take M 2 = 120 GeV, \i = 250 GeV, tan/3 
= 10, m gR3 = 150 GeV |. The sign of \x is taken so that the chargino contribution to 
the b — > sj process has the opposite sign for the SM contribution. We can see that the 

1.2 xlO -8 

1 xlO" 8 

^ 8xl(T 9 
a. 

I 6xl0" 9 
" 4xl0" 9 
2xl0 -9 

200 400 600 800 1000 1200 1400 
nWGeV) 



2x10"' 
1.5x10"' 

u 

A 

3 1x10"' 
pq 

5x10"' 



200 400 600 800 1000 1200 1400 

Figure 1: Dependence of Br{r — > /jlj) and Br(fi — > e-y) on m£ Rl2 

branching ratios become zero at the point where mg Ri is 150 GeV. This is because the 
flavor violating source (m^) is vanished, since right-handed slepton masses are univer- 
sal. Since we do not take account of the slepton LR mixing or the radiatively induced 
flavor violating source in this paper, it seems that their contributions dominate actually 
at this region. On the other hand, as we can see in Fig. 1, both branching ratios are not 
vanishing in the limit rrig Ri 3> . Let us explain the reason, using the mass insertion 
method. Both r — ► /ry and \i — * ey diagrams (a) and (b) include the first or the second 
generation right-handed slepton in the loop. However, in the limit ^e Rl2 ^ m e R3 , the 
internal lines can be estimated as 

2 This parameter set has been adopted to avoid the weak scale instability in the paper [T^ in which 
we discussed the little hierarchy problem. 
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e «3 eR l,2 e R3 



for r — > /ry, 

_ JrnUn __e R _ _ 1 _ m 2 )A 5_L_ Q ? ( 2 1) 



jUfl ( m miz) 23 Tr { m nrix) 31 2 2\2\5 ^ ^ m, ^ OO 



-X- 



(m -m 2 ) 2 A 5 ^ — (22) 



for fi — > e7. When rrie Rl2 m e R3 , the flavor changing off-diagonal entries become also 
large. As a result, both transition rates remain finite. From Fig. 1, the branching ratios 
no longer depend on the first two generation masses in me Ri 2 > 600 GeV region. 

In Figs. 2 and 4, we show the /i dependence of Brir — > ^7) and Br{^i — > cy). Here, 
we take m e - Jii = 1000 GeV and the other parameters are taken as in the previous one. 
We find that around fi ~ 20, —1000 GeV the branching ratios become very small. To 
understand this behavior, we should remind that neutralino amplitude consists of two 
parts. One is proportional to m^., the other is proportional to M^o, they correspond to 
diagram (a) and diagram (b), respectively. In Figs. 3 and 5, the dotted line represents 
the amplitude which are proportional to and the solid line represents the amplitude 
which is proportional to M%o. As we can see from these figures, the mi. part depends 
hardly on fi. This is because this part comes mainly from the diagram (a) in which 
higssinos do not contribute. On the other hand, the M^o part which comes mainly from 
the diagram (b) depends on /z. This part changes its sign from /i < to fx > and is 
decreasing as increasing /i due to the decoupling. Basically, contribution of diagram (b) 
is larger than that of diagram (a) due to tan/3 enhancement. However, when diagram (b) 
changes its sign, or it goes to decouple as increasing there are two points where the 
magnitude of diagram (b) becomes smaller than that of diagram (a). If each diagrams 
have different sign in these points, they are cancelled. Now, we could understand the 
vanishing point in Figs. 2 and 4. Note that the cancellation happens mainly in the 
fi < region. Consequently, the branching ratios do not strongly depend on [i in the 
positive /i region where the chargino contribution to b — > S7 has opposite sign for the SM 
contribution and the charged Higgs contribution to b — > sj. 

In Figs. 6, we show the M 2 dependence of the branching ratios when /i = 250 GeV. 
Here, we take the same parameters as those in the previous figures. We can see the point 
around Mi ~ GeV in which the branching ratios are vanishing. This is because of the 
cancellation between the contributions from the diagram (a) and (b) as in Figs. 2 and 4. 
The region with M 2 < and fi > corresponds to the region M 2 > and /i < 0. 

Figs. 7 and 8 show the r — > ^7 and fi — > c-f branching ratios for tan f3 and right- 
handed stau mass when M 2 = 120 GeV. We take the other parameters to be the same 
as in the previous figures. In Fig. 7, the dashed line represents Br(r — * ^7) = 4.5 x 10~ 8 
and the left side of the dashed line is excluded by current experiment. The solid lines 
represents Br{r — > /17) = 1 x 10~ 8 , 10~ 9 , 10~ 10 , 10~ n from the up-left side. It is found 
that if right-handed stau mass is about 100 GeV, tan/3 cannot be taken larger than 9. 
The region where 4.5 x 10~ 8 > Br(r — > ^7) > 5 x 10 -9 is the target region of the B- 
factory or super B-factory[T8]. If the right-handed stau is light (m^ < 200 GeV, when 
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tfGeV) 



Figure 2: n dependence of Br(r — > ^7) 




-1000 -500 500 1000 

MGeV) 



Figure 3: The solid line represents the /i dependence of the r — > ^7 neutralino amplitude 
which is proportional to m^o, and the dotted line represents the /i dependence of the r — > ^7 
neutralino amplitude which is proportional to m T . 




-1000 -500 500 1000 

MGeV) 



Figure 4: The \i dependence of Br(/j, — > e7) 




-1000 - 500 500 1000 

MGeV) 



Figure 5: The solid line represents the \i dependence of the \i — > e~/ neutralino amplitude 
which is proportional to m^o, and the dotted line represents the \x dependence of the \i — > e7 
neutralino amplitude which is proportional to m^. 
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Figure 6: M2 dependence of Br(r — > fij) and 73r(ii — ► ej) 



tan /3 ~ 10), the r — ► 7/7 process can be discovered at the B factory or super B-factory. In 
Fig. 8, the solid lines represent Br(fi — > e7) = 10~ n , 10~ 12 , 10~ 13 , 10~ 14 , 10~ 15 from the 
up-left side and the dashed line represents the current limit Br(fi — > e7) < 1.2 x 10~ n . 
The region 10~ n > Br(ji -> e7) > 10~ 14 is the target region of the MEG experiment [19j. 
Thus, we find if the right-handed stau is light {rrie R3 < 420 GeV, when tan/5 ~ 10), 
fi — > e7 may be discovered at the MEG experiment. We calculated these figures with the 
negative fi = —250 GeV. The results are similar to the Figs. 7 and 8. However, when the 
jj, is negative, the chargino contribution to the b —>■ 57 process has the same sign as the 
SM contribution and the charged Higgs contribution. Therefore, other large contribution 
to the b — >• S7 process with opposite sign is required. 

Comparing these figures, we find that the r — ► ^7 process gives more severe constraints 
to this scenario than the \i — *■ process at present. Before the MEG experiment finds 
the signal of the \l — > ej, the r — > fi^y events can be found in the B-factory. 

4 Conclusion 

In this paper, we have analyzed the lepton flavor violating r — ► /ry and fi — > processes 
in SUSY GUT model in which the third generation sfermion masses, 777,3 of 10 of SU(5) 
have different masses from the others, m. Such mass spectrum is consistent with the 
present constraints from the various FCNC processes. And it is interesting that such mass 
spectrum is naturally derived from the GUT models with E 6 x Horizontal Symmetry. 

This mass spectrum allowed us to take very large value for m without destabilizing the 
weak scale in the MSSM, because stop masses are determined by m 3 . We have considered 
a situation where m 3 is of order the weak scale while m is larger than the weak scale. In 
this case, the LFV decay rates are determined mainly by the right-handed stau mass and 
tan (3. Especially, the branching ratios of r — >• /ry and /x — >• ej strongly depend on the 
right-handed stau mass. This means that this type of soft mass matrix can be tested by 
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m gR3 (GeV) 

Figure 7: Dependence of Br(r — ► ^7) on tan/3 and m e - ffl (/i = 250 GeV). The solid lines 
represent Br(r — > ^7) = 10~ 8 , 10~ 9 , 10~ 10 , 10~ n from the up-left side to the down-right side. 
The dashed line represents the current bound. 
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m gR3 (GeV) 

Figure 8: Dependence of Br(p, — > e7) on tan/3 and m e - M (/i = 250 GeV). The solid lines 
represent Br([i — ► e7) = 10 -11 , 10~ 12 , 10 -13 , 10~ 14 , 10~ 15 from the up-left side to the down- 
right. The dashed line represents the current bound. 
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checking the relation among the LFV decay rates and the right-handed stau mass which 
may be measured in the direct search, though the ambiguity due to the 0(1) coefficients 
still exists. As we can see from Figs. 7 and 8, if the right-handed stau is lighter than 420 
GeV (when tan/5 ~ 10), fi — > ej process may be discovered at the MEG experiment. 
It may be interesting that the B-factory can find the r —>■ fij events before the MEG 
experiment find the /i — > events because the present constraints from r — > fi'-f is more 
severe than that from the fi — > process. 

Moreover, this model provides some interesting predictions. First, this model says 
that the final state lepton tends to have the right-handed chirality, because the inter- 
mediate state must have the right-handed chirality to pick up the flavor violating mass 
matrix entries. On the other hand, the LFV processes induced radiatively via large Dirac 
neutrino Yukawa interactions predicts an opposite chirality (left-handed) for the final 
state lepton. This difference can be checked by measuring the angular distribution of 
the final lepton for the initial lepton's spin. Second, this model predicts several relations 
among r — > /17, /1 — > ej, and r — ► as Br(r — > 717) ~ 10 3 x [C(l)] 4 5r(/i — > ej) ~ 
20 x [0(l)] A Br(T — > e , -f). This enhancement factor of the former relation can be under- 
stood as Brir — > evv) x A 4 /A 10 and that of the latter as A~ 2 . It is important to test 
these predictions in future experiments. 
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Appendix 

In this appendix, we give our notation and conventions adopted in Section 3, which are 
basically the same as in the papers [201 EI] ■ 

First we consider the sfermions. Below the weak scale, the mass matrix of the charged 
slepton is given in the super-CKM basis as 



is a diagonalizing matrix given in Eq. ©, and m e is the diagonal charged lepton mass 
matrix. Here we assumed that the A-term is proportional to the Yukawa matrix. We 
define the 6x6 unitary matrix T e which rotates the left(right)-handed slepton fields into 
mass eigenstate (k — 1, . . . , 6) as 




(23) 




r e (mg)r^ = (diagonal), 



e L (R) 



eL(R) e ' 



where 6x3 matrices T £ l(r) are defined as 



r 



<ki 

eL — 



r 



r 



eR ~ 



e 



(25) 
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Now we turn to the neutralinos. The mass term of the neutralino sector is given by 



/7\ 

Z 

H° d 



where 



/ Mi cos 2 9 W + M 2 sin 2 9 W cos 9 W sin 9 W (M 2 - M x ) 

;os 9 W sin W (M 2 -Mi) M x sin 2 9 W + M 2 cos 2 9 W 
mz cos (3 

\ — sin (3 



This neutralino mass matrix can be diagonalized by a real orthogonal matrix N: 



+ h.c. , 










cos P 


-m^ sin /? 













(26) 



NM^oN T = (diagonal) 
The mass eigenstates xl (i — 1> ••• ) 4) are given by 



(x°A 






xl 
xi 


= N 


z 


\x°J 




\H°J 



We use the following vertices. 



(-k 



where 



^OeL 



^OeR 

n 0eL 
jrjki 
n 0eR 



-sm9 w W u + 

-sin0 w N J ' 1 + 

Wi\T GeL Y e ) u /g, 
N*i 3 (T 0eR Y e ) k y g . 



1 



cos 9 W 2 
sin 2 ^ Ni2 



+ sin 2 9 w )N j2 



eR ) 



-pki 
1 eL ) 



Finally, we show the functions F 2 and F 4 which appear in eq. (1181) as 

1 



F 2 (x) 
FAx) 



12(2 - l) 4 

1 

2(x - I) 



(2x 3 + 3x 2 — 6x + 1 — 6x 2 logx) 
x 2 — 1 — 2x logo;) . 



/ 

(27) 

(28) 



(29) 



(32) 

(33) 

(34) 
(35) 



(36) 
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